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SUMMARY 

A multiwavelength technique has been developed and applied to measure the temperatures 
of nongray surfaces. The instruments required are a spectral radiometer, a dedicated auxiliary 
radiation source, and a computer. In general, three radiation spectra are recorded: (1) spectrum 
Sq of the auxiliary radiation source, (2) spectrum Sj of the surface-emitted radiation, and 
(3) spectrum Sjj, the sum of the radiation of Sj plus the reflected radiation due to the 
incidence of the auxiliary radiation source on the surface. Subtracting spectrum Sj from 
spectrum Sjj yields the reflection spectrum resulting from the incident radiation. From these 
spectra, a quantity z(A) is derived and is related to the reflectivity r(A) by r(A) = z(A)/f, 
where f is a constant. Spectrum Sj is represented mathematically as the product of a 
wavelength-dependent emissivity obtained from KirchhofFs law and a Planck function of tem- 
perature T. Application of two-variable (A and z), nonlinear, least-squares curve- fitting 
computer software to fit spectrum Sj to this mathematical expression yielded the surface 
temperature. This technique also measured the spectral reflectivity and emissivity of the 
surface. Instrumentation necessary to extend measurement to elevated temperatures and in the 
presence of reflective interference is discussed. 


INTRODUCTION 

Traditional pyrometry uses either a one-color or a two-color method. In one-color pyro- 
metry, temperature determination requires knowing the emissivity of the surface. In the range 
of its design operating temperatures (usually in the neighborhood of 1000 K), the fractional error 
in temperature measurement due to error in emissivity is only about one-tenth of that due to 
error in pyrometer radiation signals (ref. 1). Nevertheless, accurate temperature measurement 
would still require more accurate values of emissivity. Accurate temperature measurement with- 
out prior knowledge of emissivity is possible if a functional dependence of emissivity with 
wavelength is assumed (refs. 1 and 2). In two-color pyrometry, more accurate temperature mea- 
surement is possible without exact knowledge of emissivity by assuming that the emissivity is 
the same at the two wavelength bands where measurements are made (ref. 3). However, neither 
method is capable of correcting for the interference of reflected radiation. 

At NASA Lewis, multiwavelength pyrometry has been used to measure surface 
temperature and to correct for the interference radiation from an extraneous source (refs. 4 
and 5). The measured pyrometry radiation spectrum in general covers a broad wavelength 
region. In these previous studies the surface under consideration was assumed to be gray 
(emissivity independent of wavelength). In actual application, many bodies such as ceramic 
materials are nongray. For these materials, the emissivity is known to be wavelength and 
temperature dependent. It is small at short wavelengths and increases to almost unity at 
wavelengths around 8 or 10 /im, and then decreases and increases in some manner at longer 
wavelengths (ref. 6). Thus, conventional pyrometry optimized only for black or gray bodies 



would produce unavoidable error when applied to nongray bodies (ref. 7). Even for a surface 
with known emissivity, the value may still vary over time. Accurate temperature measurement 
requires measuring the emissivity simultaneously. 

As described in this report, multiwavelength curve-fitting pyrometry can be extended to 
measure accurately the surface temperature of a nongray surface by simultaneously measuring 
the emissivity. It can also accommodate the more general and important case in which 
extraneous interfering radiation is present. This will be described in another report (ref. 8). 
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THEORY 

The basic components of this pyrometer are a spectral radiometer, an auxiliary radiation 
source which can be turned on and off (either manually or by mechanical chopping), and a 
computer for data acquisition, data processing, and numerical calculation. 

Three spectra are required: 

(1) S 0 (A), the spectrum of the auxiliary radiation source (can be accumulated prior to real time 
measurement) 

(2) Sj(A), the spect rum of the surface radiation 

(3) S n (A), the spectrum resulting from the superposition of radiation emitted from the surface 
and radiation from the auxiliary radiation source reflected by the surface 

Spectrum S 0 (A) is the direct unobstructed spectrum of the auxiliary radiation source in 
the wavelength region of interest obtained in the geometry shown in figure 1 using the spectral 
radiometer of a multiwavelength pyrometer. 

Spectrum Sj(A) is obtained in a typical pyrometer geometry shown in figure 2 with the 
auxiliary radiation source turned off. This is the emission spectrum of the nongray surface 
described by a Planck function of temperature T modified by the wavelength-dependent 
emissivity of the nongray surface. 

Spectrum S n (A) is also obtained in the geometry shown in Figure 2 but with the auxiliary 
radiation source turned on. Thus S n (A) is the spectrum Sj(A) plus reflected radiation due to 
the auxiliary radiation source. Subtracting spectrum Sj(A) from S n (A) produces the reflected 
spectrum S m (A) = S n (A) - Sj(A) due to the incidence of the auxiliary radiation source on the 
measured surface. 


The measured reflectivity at any wavelength is the ratio 

... S, U (A) S U (A)-S,(A) 

zf A) = = 

S 0 (A) S 0 (A) 
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Because S Q (A) is measured on a direct path in figure 1, and Sjjj(A) is measured on a reflected 
path in figure 2, their ratio z(A) includes a constant factor f due to different conditions along 
the two paths. Two conditions are the geometry of the optical beams and the nonspecular 
reflecting surface. Thus, the true reflectivity is 


r(A) = 

f 


( 2 ) 


where f will be determined. 

KirchhofFs law is applied to yield 


e(A, /9, 0 , T e ) = 1 - r(A, 0, T e ) 


in which emissivity and reflectivity generally possess spectral, angular, and temperature depen- 
dencies. The symbols /? and 6 refer to the polar angular coordinate angles specifying a 
direction with reference to a suitably chosen coordinate system, and T e is the temperature of 
the reflecting surface. In pyrometry application, the angular dependence disappears because the 
optics of a pyrometer’s detector selects signals in a narrowly defined direction. As a result, only 
the wavelength dependence needs to be considered at a particular temperature. The emissivity 
of the surface is therefore given by 


e(A) = 1 - r(A) 

= 1 _ ^) 
f 


( 3 ) 


In principle, the constant f can be determined from the geometry of the experimental setup. 
However, it is more conveniently determined by curve fitting. 

The mathematical expression for the surface emission radiation spectrum Sj(A) is a 
Planck function of temperature T modified by the wavelength-dependent emissivity so that 

C, , 

_ ( 4 ) 
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Substituting equation (3) gives 


Si(A) = 1 




f ) 


exp 


AT 


( 5 ) 


where c l and c 2 are the first and second radiation constants, and f and T are unknown 
parameters. The quantity z is obtained using the term-by-term spectrum subtraction and 
division procedures described previously in this section. The adjustable parameters in the 
expression, f and T, are determined by applying two-variable (A and z) least-squares curve 
fitting to fit the spectrum Sj(A) according to equation (5). 

RESULTS 

Multiwavelength pyrometry was used to measure the surface temperature of a silicon 
carbide (SiC) ceramic. Silicon carbide is a common ceramic material with well-documented 
wavelength-dependent emissivity and reflectivity (refs. 6 and 9). The surface of a silicon carbide 
wafer sample measuring 3 by 25 by 50 mm (0.125 by 1 by 2 in.) was nominally polished to a 
10-pm finish. The temperature of the sample was raised above ambient to about 80 °C by 
placing it in front of a blackbody furnace and allowing to equilibrate. The blackbody furnace 
temperature was regulated by a temperature controller with ±0.5 °C long-term stability. The 
emission spectrum was obtained using a spectral radiometer in the wavelength region 2.5 to 
14.5 pm. 

The spectral radiometer that was used produced a spectrum either in voltages or in radia- 
tion energy units. When the spectrum was in volts, the value was a function of wavelength 
depending on the intensity of the radiation it detected and the detector responsivity. The 
voltage spectrum was converted into a radiation energy spectrum when internally divided by the 
spectrometer’s response function which was obtained during routine calibrations performed 
periodically. 

The emission spectrum (in radiation energy units) is shown in figure 3 as a function of 
wavelength. The deviation of the emission spectrum from a blackbody Planck curve correspond- 
ing to the SiC surface temperature is clearly visible. This deviation is due to the wavelength- 
dependent emissivity. The spikes in the wavelength region 6 to 7 pm are artifacts due to 
absorption of water and carbon dioxide in the atmosphere. However, the dip at 12 pm is due to 
a sharp rise (decrease) in reflectivity (emissivity). This increase in reflectivity is explained 
by the classical dispersion theory of crystals. Since the radiation energy spectrum is a secondary 
quantity derived from the primary voltage spectrum, the primary voltage data are used in the 
following analysis which describes how the reflectivity is obtained. 
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Figure 4 shows the surface emission spectrum Sj of the SiC surface (in volts) obtained by 
the spectrometer according to the geometry shown in figure 2. Also shown is the spectrum Sn 
(also in volts) obtained when the surface is illuminated by an auxiliary source, a 20-W com- 
mercially available infrared lamp. The lamp was regulated by a constant-current power supply 
such that its typical light output RMS ripple was 0.05 percent. The difference between these 
two curves gives the reflected spectrum (in volts) due to the auxiliary source incident on the SiC 
surface (fig. 5). The spectrum S Q of the unobstructed auxiliary light source acquired with the 
geometry shown in figure 1 is shown in figure 6. The value z is calculated according to equa- 
tion (1) and shown in figure 7. This is an experimentally derived quantity related to the emis- 
sivity. The quantity z(A) is obtained by using the term-by-term spectrum subtraction and 
division procedures described previously. The adjustable parameters (f and T) in the expres- 
sion are determined by the application of two-variable (A and z), nonlinear, least-squares curve 
fitting to fit the data of spectrum Sj(A) to equation (5). This curve-fitting software is commer- 
cially available. 

The temperature determined by curve fitting was 363 K. The good agreement of the 
fitted curve with the experimental spectrum is shown in figure 8. A 0.125-mm (5-mil) type K 
(Chromel-Alumel) thermocouple was used such that its junction made as good a contact with 
the SiC surface as possible. It measured a surface temperature of 360.2 K, in very good 
agreement with the multiwavelength pyrometer determined value. 

The least-squares determined value of f was 0.0248. It was used to calculate the reflec- 
tivity according to equation (2) and the emissivity according to equation (3). These results are 
shown in figures 9 and 10, respectively. The reflectivity’s spectral dependence agrees qualita- 
tively with published reflectivity from a polished, optic-axis-oriented crystalline SiC sample using 
polarized radiations in the 1- to 20-jim wavelength region (ref. 9). The reflectivity r(A) is satis- 
factorily described by the classical dispersion theory of crystals given by 

r(A) = ~ ^ + kZ (8) 

(n + l) 2 + k 2 


with 


n 2 = (1/2) 


U* + 4 


( \2 
a 


11/2 


+ f 


( 9 ) 
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k 2 = (1/2) 


e 2 + 41 
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f = <0 + **X 


( 11 ) 
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characterized by c Q , i/ 0 , p y and 7, which are the high-frequency dielectric constant, the 
resonance frequency, the resonance width, and the strength of the resonance, respectively. These 
are also related to the susceptibility conductivity a, index of refraction n, and extinction 
coefficient k. For ease of calculation, v is the measured frequency divided by i/ Q , and a is the 
conductivity divided by i/ Q . The reflectivity so calculated using values of z/ Q = 2.75xl0 13 Hz, 
p — 0.029, 7 = 0.15, and e Q = 2.4 is shown in figure 9 together with the experimentally 
measured results. In this way the experimentally determined reflectivity data is fitted to 
equation (8). 

At temperatures well above ambient, the difference between the signals of the radiation 
spectra with and without the auxiliary source is not easy to measure. However, the signal 
corresponding to that difference can be readily measured by chopping the auxiliary signal and 
using phase-sensitive detection schemes. Most spectral radiometers have external signal 
chopping capability. In this way, the spectrum of the reflected signal is readily obtained, and 
direct two-variable least-squares curve fitting will determine the temperature of the nongray 
surface. 


CONCLUSION 

Multiwavelength pyrometry can be used to determine the surface temperature of a 
nongray (emissivity depending on wavelength) surface by acquiring and analyzing three spectra, 
two of them from a nongray surface. They are the spectrum of an unobstructed auxiliary 
radiation source, the surface emission spectrum of the nongray surface, and the spectrum 
consisting of radiation due to surface emission and reflection due to the incidence of the auxiliary 
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radiation source on the surface. The reflectivity (hence emissivity) of the surface is obtained 
simultaneously during least-squares curve fitting. The reflectivity is analyzed satisfactorily using 
the classical dispersion theory of crystals. This pyrometry was demonstrated by measuring the 
surface temperature of silicon carbide, a material exhibiting wavelength-dependent emissivity. 
Extension to higher temperatures is possible by a chopping and phase-sensitive detection 
method. 
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Figure 1.— Experimental arrangement to determine 
unobstructed spectrum Sq of auxiliary radiation 
source. 
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Figure 3. — Surface emission spectrum (in arbitrary energy 
units) of SIC surface compared with Planck curve for 
same temperature. 


8X10- 3 



Wavelength, ^m 


Figure 5.— Reflected spectrum S 1H from SIC surface obtained 
by subtracting Sj from S n in figure 4. 
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Figure 2.— Experimental arrangement of pyrometry 
components. 
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Figure 4. — Surface emission spectra (In volts) of SIC surface 
(Sj, auxiliary source off; S I|t auxiliary source on). 
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Figure 6.— Unobstructed spectrum Sq (In volts) of auxiliary 
radiation source. 


8 




Measured reflectivity 



9 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information I* estimated to average 1 hour per response, Indudng the time for reviewii^^ 
gathering and maintaining the data needed, and completing and reviewing the collection of Information Send comment* [regarding this burden M«matec > 
collection of information, Including suggestions for reducing this burden, to Washington Headquarters 90W3 
Davis Highway Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Waahtngton, D C 20503. 


1. AGENCY USE ONLY (Leave blanK) 


REPORT DATE 

February 1992 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


4. TITLE AND SUBTITLE 

5. FUNDING NUMBERS 

Multiwavelength Pyrometry for Nongray Bodies 

WU-510-01 -50 

6. AUTHOR (S) 


Daniel Ng 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

E-6616 

9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESSES) 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

National Aeronautics and Space Administration 
Washington, D.C 20546-0001 

NASA TM- 105285 

11. SUPPLEMENTARY NOTES 


Responsible person, Daniel Ng, (216) 433 - 3638. 


12a. DISTRIBUTION/AVAILABILFTY STATEMENT 

12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 35 



13. ABSTRACT (Maximum 200 words) 

A multiwavelength technique has been developed and applied to measure the temperatures of nongray surfaces. The 
instruments required are a spectral radiometer, a dedicated auxiliary radiation source, and a computer. In general, 
three radiation spectra are recorded: (1) spectrum S 0 of the auxiliary radiation source, (2) spectrum Sj of the surface- 
emitted radiation, and (3) spectrum S n , the sum of the radiation of Sj plus the reflected radiation due to the incidence 
of the auxiliary radiation source on the surface. Subtracting spectrum Sj from spectrum Sjj yields the reflection 
spectrum resulting from the incident radiation. From these spectra, a quantity z(X) is derived and is related to the 
reflectivity r(X) by r(X) - z(A,)/f, where f is a constant. Spectrum Sj is represented mathematically as the product of a 
wavelength-dependent emissivity obtained from Kirchhoff s law and a Planck function of temperature T. Application 
of two- variable (X and z), nonlinear, least-squares curve-fitting computer software to fit spectrum Sj to this mathemati- 
cal expression yielded the surface temperature. This technique also measured the spectral reflectivity and emissivity 
of the surface. Instrumentation necessary to extend measurement to elevated temperatures and in the presence of 
reflective interference is discussed. 


14. SUBJECT TERMS 

Pyrometry; Curve fitting; Nongray emissivity; Reflectivity; Temperature 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURFTY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 

10 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01 -280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 














Nil, | I I i UUIIII IttU jlllMMWitiiMl u il I Mill II 






